Sub-diffusion and localization in the one dimensional trap model 
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We study a one dimensional generalization of the exponential trap model using both numerical 
simulations and analytical approximations. We obtain the asymptotic shape of the average diffusion 
front in the sub-diffusive phase. Our central result concerns the localization properties. We find 
the dynamical participation ratios to be finite, but different from their equilibrium counterparts. 
Therefore, the idea of a partial equilibrium within the limited region of space explored by the walk is 
not exact, even for long times where each site is visited a very large number of times. We discuss the 
physical origin of this discrepancy, and characterize the full distribution of dynamical weights. We 
also study two different two-time correlation functions, which exhibit different aging properties: one 
is 'sub-aging' whereas the other one shows 'full aging'; therefore two diverging time scales appear in 
this model. We give intuitive arguments and simple analytical approximations that account for these 
differences, and obtain new predictions for the asymptotic (short time and long time) behaviour of 
the scaling functions. Finally, we discuss the issue of multiple time scalings in this model. 



PACS numbers: 75.fO.Nr, 05.20.-y, 02.50.-r 
I. INTRODUCTION 

A lot of efforts have been devoted to the theoretical 
study of aging phenomena in the past decades [1-3]. Spin 
glass models, which exhibit a very rich phenomenology, 
have been widely studied theoretically both using analyt- 
ical techniques for the mean field models, or by numerical 
simulations in the finite dimensional cases. Besides these 
microscopic spin models, a simpler but phenomenologi- 
cal picture, the 'trap model', has been proposed in order 
to describe the phase space dynamics in a coarse-grained 
manner [4]. This model seems to capture, at least qual- 
itatively, some of the physics involved in the aging dy- 
namics of several systems beyond spin glasses, such as 
fragile glasses [5-7], soft glassy materials [8,9], granular 
materials [10], pinning of extended defects (such as do- 
main walls, vortices, etc.) [11]. This trap model has been 
studied mainly in its fully connected (or 'mean field') ver- 
sion [12-14], which has recently been shown to describe 
exactly the long time dynamics of the Random Energy 
Model when the distribution of trap depth is exponen- 
tial [15]. This version of the mean field model already 
exhibits a number of interesting features, such as a tran- 
sition between a stationary, 'liquid' phase, and an aging 
'glassy' phase, violation of the Fluctuation Dissipation 
Relation [16], and dynamical ultrametricity [17,18]. In 
the glassy phase, the dynamics is strongly intermittent, 
since most of the time nothing happens, whereas the ac- 
tive periods appear in bursts which become less and less 
frequent as time elapses. Several recent experiments sug- 
gest that such an intermittency is indeed present in glassy 
systems [19-21], or in atomic physics [22-24]. 

The finite dimensional generalization of this model has 
already been studied many years ago [25-28], but only 
one time quantities (not well suited to study aging) were 
considered. These aging properties were addressed only 



recently in [14,29], and, from a more rigorous point of 
view, in [30,31]. One expects that in dimensions d > 2, 
the trap model will have properties qualitatively similar 
to the fully connected case, since each site is visited by 
the walk a finite number of times. In lower dimensions 
d < 2, the correlations induced by the multiple visits of 
the walks to a given site is expected to lead to qualita- 
tive changes. It was for example shown in [29] that some 
quantities exhibit sub-aging properties, i.e. decay on a 
time scale that scale with the waiting time as with 
1/ < 1. Because of the limited number of accessible sites, 
one might also expect interesting properties such as dy- 
namical localization, which means that there is a finite 
probability that k independent particles sit on the very 
same site, even after a very long waiting time i^,. Such 
a dynamical localization was first established by Golosov 
in the context of the Sinai model [32] and extended to the 
biased case in [33], and more recently proven rigorously 
for the one dimensional trap model in [30]. 

In this paper, we present a detailed study of the one 
dimensional (non biased) trap model, using both numer- 
ical simulations and analytical approximations. In the 
first section, we focus on the scaling form of the aver- 
age "diffusion front" {p(x,t)) in the sub-diffusive, non 
Gaussian phase, for which no analytical results are (to 
our knowledge) available. We present some scaling ar- 
guments and approximation schemes to account for our 
numerical data. We then discuss the idea of partial equi- 
librium in this model, which can be explored in details 
through the distribution of dynamical weights. The mo- 
ments of this distribution are the usual 'participation ra- 
tios' that characterize the localization properties of the 
measure. Perhaps surprisingly, these localization indica- 
tors are indeed finite (as first shown in [30]), but different 
from their static counterparts. We discuss in detail the 
origin of this difference, and try to characterize quan- 
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titatively the distribution of dynamical weights. In the 
last section, we study the aging behaviour of two differ- 
ent correlation functions, which exhibit different scaling 
properties, meaning that two different time scales, and 
tw, appear in this model. We again develop intuitive ar- 
guments and simple analytical approximations to under- 
stand these differences, and obtain new predictions for 
the asymptotic behaviour of the scaling functions, which 
are found to be in excellent agreement with the numerics. 
Finally, we discuss the possible existence of multiple time 
scalings in this model (as can happen in spin glasses [17] 
or in a generalization of the trap model [29]). 

II. THE ONE DIMENSIONAL TRAP MODEL. 
ANOMALOUS DIFFUSION 

A. Definition of the model 

Consider a one dimensional lattice, and define on each 
site i a quenched random variable Ei > chosen from a 
distribution p{E). Ei has to be interpreted as the energy 
barrier that the particle (the walker) has to overcome in 
order leave the site. The dynamics is chosen to be acti- 
vated with temperature T, which means that the escape 
rate Wi of site i is given by Wi = Foe"^'/^, where Fq 
is a microscopic frequency scale. Once that particle has 
escaped the trap, it chooses one of the two neighbouring 
sites, with probability g_ for the left one and q+ = 1 — q- 
for the right one. The 'directed' case g+ = 1 is quite sim- 
ple to analyse analytically, since each since is visited once 
- see [34,28,33]. The case q+ = 1/2 that we study in the 
following is much more subtle since each site is visited a 
large number of times, inducing long range correlations 
in the hopping rates seen by the walker. (Note that as 
soon as q+ ^ 1/2, one expects the large time properties 
of the walk to be the same as in the fully directed case 
[33].) 

For the purpose of heuristic arguments and Monte- 
Carlo studies, it is interesting to study the trapping 
time T of the particle on each site. Once the transi- 
tion rates Wi arc given, r is a random variable with a 
(site dependent) distribution ^(t) = Wi e~'^''^, of mean 
Ti = . If we choose an exponential density of trap 
depths, p(E) — ^e^^/'^n, then the distribution of r^'s 
over the different sites is a power law: 

where /i — TjTg is the reduced temperature, and tq = 
Fq ^. For T > Tg, this distribution has a finite average 
value (r) = fo/(/i — 1). This corresponds to usual diffu- 
sion and stationary dynamics, with a diffusion constant 
D ~ /{t), where a is the lattice spacing. On the con- 
trary, for T < Tg, the first moment of the distribution 



diverges, diffusion becomes anomaloiis and aging effects 
are expected. A dynamical phase transition takes place 
at Tg, as in the fully connected model. However, new 
properties emerging from the non trivial spatial struc- 
ture of the model are expected. 

B. Disorder induced sub-diffusion: A scaling 
argument 

We first give a simple scaling argument (proposed in 
[25,35,28]) that yields a sub-diffusive behaviour for the 
one dimensional trap model introduced above. In the 
following, we shall take a is the unit of length, as well as 
Fq ^ as the time unit. Roughly speaking, a typical ran- 
dom walk starting from a given initial site has visited, 
after N steps, of the order of VN sites (which implies 
a typical displacement ^ ~ VN)- Each site is visited 
around times. So the time t elapsed can be written 
as: 

Since the sum of M independent random variables dis- 
tributed according to Eq. (1) grows as M>^, we get: 

t^^/N''^^ r^C^i. (3) 

Inverting this relation leads to the following sub-diffusive 
behaviour: 

C(t)~tW. (4) 

This result was also obtained by Machta [27], using Real 
Space Renormalization Group arguments. The same be- 
haviour also holds for the random harrier model with a 
broad distribution of barrier heights which, in one dimen- 
sion, is expected to be equivalent to the trap model, as 
far as diffusion properties are concerned [28]. For this 
model, the average probability of being on the initial site 
can be exactly computed, and decays as 1/^1+'' = 1/^(^)5 
in agreement with the above result. The exponent is 
also in very good agreement with numerical results [29] . 

The case = 1 is special since logarithmic corrections 
come into play. Extending the above argument leads to: 

^^)-\[^t (/^ = 1), (5) 

whereas for > 1 one recovers £^{t) ~ ^/t. 

Calling p{x = ia,t) = Pi {t) the probability to be at a 
distance x from the starting point after time t, one ex- 
pects the disordered average diffusion front {p{x, t))r to 
take for large times the following scaling form: 
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where ^{t) is given by Eq. (4), and (. . .),- stands for the 
average over the quenched trapping times r^. However, 
the full scaling function /(.) is, to our knowledge, not 
known. Only the value of /(O) in the dual 'barrier' model 
was obtained in [26] . Before studying more subtle issues, 
we have investigated this question both analytically and 
numerically. 



C. The average diffusion front 

For a system without disorder, or in the case /x > 1 
where the average trapping time (r) is finite, the Central 
Limit Theorem tells us immediately that the diffusion 
front becomes Gaussian at large times: 



{p{x,t))r 



(7) 



where D = / {t) is the diffusion constant, a being the 
lattice spacing, so that D (x (/U — 1) when /x — > 1+. In the 
case /X < 1, a modified space-time scaling is expected, 
as argued in the previous subsection, as well as a non 
Gaussian diffusion front. 




FIG. 1. Plot of i{t){p{x,t)) versus x/^{t) for different 
temperatures and different times. Data were obtained by 
Monte-Carlo simulations on the model with fixed trapping 
times. Upper curves: = 0.2, and t = 10^ (o), 10* (A), 
10^° (+); middle curves: fj. = 0.5 and t = 10^ (□), 10* (<), 



10 
10 



^ ( X ) ; lower curves: 
(V). 



H = 0.9 and t = (>), 10* (O), 



We have developed simple approximation schemes 
(that we expect to become exact in the limits fi 1 
and t oo) to compute {p(x,t))r analytically. The 
calculations are reported in Appendix A. We find that 
{p{x, t))r can indeed be written as Eq. (6) with £^{t) given 
by Eq. (4). The asymptotic shape of /(C = x/^{t)) can 
furthermore be computed in the limits C ^ and C ^ oo. 
We find: 



/(C)«/oo|Crexp(-6|C|^) Id -00 (8) 

/(c)«/o-/iicr-/2icr Ki-o, (9) 



where a = (/z — l)/2, /3 = 1 + /U and 7 = min(2, 1 + 2/z). 
The constants /o, /i, /2, /oo and b are n dependent num- 
bers that we can also compute in the Appendix (see 

Eqs (83) and (90)). 

In the case /x = 1, we find, using the same approxima- 
tion (which we now believe is exact), that C is given by 

Xy/^^cjt w x^ \ \nt/t, and that /(C) is exactly Gaus- 
sian, as for the 'normal' case /i > 1. More precisely, one 
finds for the diffusion front, in the limit t 00: 



{p{x,t))r ^ 




exp 



X , 

\nt 

At 



(10) 



We have tested numerically the validity of the scaling re- 
lation (6), for several values of /u. The plot of ^(f) (j){x, t)) 
as a function of x/^{t), for different values of t shows a 
rather good collapse. The curves collapse well for /i = 0.2 
and 0.5, even if for /x = 0.2 data is more noisy. However, 
finite time corrections become stronger as /i approaches 
1. This is expected: sub-leading corrections to scaling 
can be shown to become negligible only in the limit where 
f 1+M ^ 1. For ji = 0.9 and t = 10 , however, this param- 
eter is only « 1.8. Therefore, we expect that the fi = 0.9 
data will actually be strongly affected by the vicinity of 
/i = 1, which plays the role of a critical point. 




FIG. 2. Plot of |Cp"'''''^fl(C) versus 1(1^+" for /tt = 0.5 (O), 
0.8 (o) and 0.9 (A), obtained by an infinite time extrapolation 
of the Monte-Carlo data. The analytical prediction, using the 
approximation valid for fi close to 1, is also shown for the same 
values of yu (lines). The predicted exponent ICI^^*^ is in good 
agreement with the numerics, since data appear to be linear 
in this representation. The prediction for b - sec Eq. (8) - is 
in best agreement with the Monte-Carlo data for fi — 0.8. 
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In Fig. 2 wc show the scaling functions for /i = 0.5, 0.8 
and 0.9 obtained by extrapolating to t = oo the scaling 
curves obtained at finite t. We actually plot /(C)/C" as a 
function of in a semi- log plot, in order to test directly 
the asymptotic form given by Eq. (8). Note that the ap- 
proximation is supposed to be valid only for /x close to 
1 , but seems to work well even for rather small values of 
12, like 12 = 0.5. Also shown are the analytic predictions, 
with the computed numerical values of the constants foo 
and b. We see that the agreement is quite reasonable, and 
actually suggests that the value /? = 1 -|- /x is probably 
exact. For fj, = 0.9, critical corrections become important 
and the predicted slope is not as good as for /x = 0.8, but 
the exponent I + ii seems to be correct. 

Data corresponding to = 1 are shown in Fig. 3. It 
is actually necessary to take into account a finite time 
correction in this case, replacing Int by In(ri), where T 
is an unknown constant which has to be fitted. This cor- 
rection is natural, as can be shown by the structure of 
the sub- leading terms. Fitting F on data corresponding 
to < = 10'* leads to F ~ 1.64. One sees that for large x, 
the slope is then found to be the same for the three times 
simulated. 





t=io' 


• 


1=10" 


A 


t=10' 




Prediction 




FIG. 3. Plot of ^c{t){p{x,t))r versus [x/^c{t)f for n = 1 
and times t = 10^ (o), 10* (•) and W'^ (h); ^a{t) is the criti- 
cal coherence length, defined as ^dt) = [t/\n{Tt)]^ , where F 
is fitted on the t = 10* curve. The data points seem to con- 
verge at large times towards the asymptotic Gaussian form 
predicted by Eq. (10) (line). 



III. 'PARTIAL EQUILIBRIUM' AND 
LOCALIZATION 



The one dimensional diffusion problem is interesting 
because, as mentioned above, each site is visited by the 
walk a large number of times. A natural idea is therefore 
that at time t, the probability Pi{t) to find the particle 
at site i should be very similar to the equilibrium distri- 
bution restricted to an interval of finite length oc ^(t). 
More precisely, we can expect that Pi{t) can be written 
on the following 'quasi-equilibrium' form P?^{t): 



p^{t)^pr{t) 



(11) 



J2 9^it)' 



.Ei/T 



where the 'form factors' gi{t) are slowly varying and de- 
cay on the scale of ^(t). (Note that the energy barrier 
Ei > Q is the opposite of the energy of the site.) This 
idea of 'partial equilibrium' is actually quite general and 
is often advocated in the context of glassy dynamics: al- 
though the system is o\it of equilibrium, one may think 
of its state at time t as of a partial equilibrium restricted 
to the region of phase space that it has explored up to 
time t, see e.g. [36-38]. This idea of partial equilibrium 
was introduced and used quantitatively in the context of 
random walk models in [29] . 

In this section, we want to discuss this issue in some 
details. It turns out that the full statistics of Pf^ can 
be worked out in the limit where ^{t) oo, and can 
be compared to the corresponding statistics of Pj (t) that 
we determine numerically. Perhaps surprisingly, we find 
that these statistics differ significantly even in the long 
time limit, meaning that the out of equilibrium problem 
never approaches a 'quasi-equilibrium' regime. 



A. Participation ratios and localization 

In order to investigate the statistical property of a ran- 
dom probability measure (such as Pi{t) or one can 
introduce the following distribution; 

ip{P) = {Y,PS{P-Pi))r 0<P<1 (12) 



We have also tested the small ( region. For /x = 0.5, 
the vICI singularity predicted by our approximation - see 
Eq. (8) - is rather convincing. However, as /i increases 
towards 1, the coefficient /i of jCj'* decreases towards 
zero. The next leading term becomes important and one 
indeed observes an effective singularity with an exponent 
intermediate between fi and 2: we find this exponent to 
be w 1.6 for fi = 0.8 and w 1.8 for = 0.9. 



which is defined in such a way as to give a small weight 
to the very large number of sites with small energies, 
in order to evidence the statistics of the deeper traps 
present in the system. This distribution ip{P) is normal- 
ized, since: 



[\{P)dP={Y,Pi)r = l 

Jo , 



(13) 
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The moments of this distribution are related to the so- 
called inverse participation ratios Yk'. 



Dynamical localization and weak ergodicity 
breaking 



[' P''-^^{P)dP = {Y,Pt)r=Yk 

Jo . 



(14) 



These participation ratios have an interesting interpreta- 
tion: if Yh remains finite for A; > 1 as the number of terms 
in the sum diverges, one speaks of localization, since a 
finite fraction of the number of particles remain concen- 
trated on a finite number of sites, even in the limit of an 
infinite number of available sites. The Yk where intro- 
duced in the context of electronic localization [39] and in 
spin glass theory [40] , and studied in several other prob- 
lems [41,42]. Note that in the limit k ^ 1, {Yk-l)/{l-k) 
becomes the statistical entropy of the measure Pi. 

In equilibrium, and for integer values of k, Yk can be 
interpreted as follows: suppose one chooses at random k 
particles with their corresponding equilibrium weight, Yk 
is the probability to find them all at the same site. Cor- 
respondingly, for the out of equilibrium situation, Yk is 
the probability to find k particles (that all started at the 
same site) clustered together on the same site at time 
t. Obviously Yk can only be non zero if some effective 
attraction exists between the particles. In the case of 
disordered systems, this attraction is induced by the dis- 
ordered environment, where (non interacting) particles 
condense into particularly favorable sites. 

For the problem at hand, the quasi-cquilibrium value 
Yj5* of the participation ratios can be computed, using 
for instance auxiliary integrals [41] (see also [43,44]). In- 
terestingly, the detailed shape of the 'form factors' c;,; in 
Eq. (11) does not matter, and the quasi-equilibrium re- 
sults yj?^ coincide with the equilibrium values Y^'^ . For 
/i > 1, Y^"^ tends to when £_{t) goes to oc', whereas for 
/Lt < 1 it converges for large ^(t) to a finite value: 



T{k - ^i) 



identical to that found in the Random Energy Model. 

This means that in the low temperature phase T < Tg, 
the equilibrium measure localizes over a finite set of sites. 
The corresponding equilibrium distribution of weights is 
given by: 

^-''(^) = r(i-!.)r(,) ^"^(^-^)^"^ 

The precise question we want to ask, in order to test 
the partial equilibrium idea, is whether or not the dy- 
namical Yk{t) approach, in the long time limit, the equi- 
librium Y^'^, computed for a large system size L. This 
can be also seen as a question about the commutation 
of the two limits t ^ oo and L ^ oo (see e.g. [3]). The 
equilibrium case corresponds to taking t ^ oo first, at 
fixed L, and then taking L to infinity. The out of equi- 
librium case, on the other hand, corresponds to taking 
L = oo from the outset and let t — > oo. 



Let us now turn to the dynamical localization prop- 
erties, starting from a localized initial condition, Pi{t = 
0) = Sifi. It has been shown by Fontes et al [30] that 
the random walk process with a diverging local mean 
trapping time converges, up to a space-time rescaling, 
to a stationary process. Consequently all spatially in- 
tegrated (one time) quantities like participation ratios 
converge to asymptotic values at large time, which are 
a priori different from the equilibrium ones. Unfortu- 
nately, this mathematical approach has not been able 
yet to predict the corresponding numerical values. We 
have computed numerically, using a simple Monte-Carlo 
method, the time dependence of Yk{t) for several values 
of fc (fc = 2, |,3, |,4); see Appendix B for technical de- 
tails. Our simulations confirm the convergence of Yk{t) 
towards a limiting value for large t, and show that these 
asymptotic values are indeed different from the equilib- 
rium ones. 




FIG. 4. Test of the convergence ofY^it) towards an out of 
equilibrium value, for jj, = 0.5; Y'zit) is plotted for different 
small sizes L = 2N + 1 of the system, so tliat equilibration 
can be reached within simulation time. One can see the onset 
of a plateau at a value lower that Y^'' = 0.5 (the equilibrium 
value for L — » ooj. The inset shows that the curves collapse 
if time is rescaled by the equilibration time terg oc N^^'^'-^^^^, 
so that the plateau corresponds to a true out of equilibrium 
effect, and not to an initial transient. 



In order to evidence the convergence of Yk (t) towards 
different asymptotic values depending on the order of lim- 
its t ^ oo and L ^ oo, we have first studied small sys- 
tems for different sizes L = 2N -|- 1, in the case fx = 0.5. 
Fig. 4 shows the onset of a clear plateau at a value 



5 



y^^"(L) smaller than the value predicted by Eq. (15), 
¥2"^ = 0.5 (for L 00), before a cross-over towards the 
equilibrium regime. Rescaling the time coordinate by a 
factor A^(i+^)/^ (corresponding to the equilibration time 
terg of the System), the data collapse rather well, at least 
in the cross-over region. This shows that the plateau in- 
deed corresponds to the onset of an out of equilibrium 
steady state regime, when the diffusion length is much 
smaller than the size of the system. The cross-over ap- 
pears when the two lengths become comparable. 

In order to study the asymptotic value y^'^^" = 
Y2^^{L 00), we have simulated systems of very large 
sizes L. However, the temporal convergence of Y2{t) is 
very slow and some infinite time extrapolation procedure 
is needed. As illustrated in the inset of Fig. 5 for fc = 2 
and fi = 0.5, we have assumed a power-law convergence 
of Yfc , of the form Yk (t) = y/^" + Air'' with 3 fitting pa- 
rameters Yl^^'\ A and a, which was found to work rather 
well. However, for fj, close to 1, the fitting parameter 
Y^^^ becomes very sensitive to the choice of the time 
interval used to fit the data, and as a result, error bars 
become larger. 




FIG. 5. Plot of Y2 (upper symbols) and Y^ ^" (lower 
symbols) vs. ji. The equilibrium functions Y2''{ij,) (full line) 
and Y^''{ij,) (dashed line) are shown for comparison. One 

clearly sees that localization is weaker than in the equilib- 
rium situation. In particular, the participation ratios seem 
to converge to a zero temperature limit which is less than 1. 
The stars and the dotted lines correspond to the prediction 
of a simple model given below, Eq. (39). Inset: fit ofY2{t) 
using the functional form Y2{t) = y/"" + At''', for fi = 0.5. 
Only 1 Monte-Carlo point out of 12 is shown, for clarity. 

Fig. 5 shows the extrapolated ¥2^'^ and Y^'"'^ as a 



function of fi, and compares it to the equilibrium rela- 
tion Y2''{n) = l-H, and Y^^fi) = (1 - At)(l - m/2). It 
appears that the dynamical localization is weaker than 
in equilibrium. In particular, y^''^" and Y^^" converge 
to a value smaller than 1 when /x goes to 0. We shall 
argue below that ^^^"(/i = 0) = 2/(fc + 1), whereas 
Y^'^ip = 0) = 1. In the other limit, /i ^ 1, it will also 
be argued in the next section that Y^^"' vanishes linearly 
with /It. This is indeed compatible with the numerical 
data, although other functional dependence might also 
be compatible since the error bars are large in this range 
of /X. We have also shown in Fig. 5 the prediction of a 
simple argument given in section IV D below, which sug- 
gests Y^^"{fx) = 2Yjf'(/x)/(fc-|-l), which is in rather good 
agreement with the numerical results. 

Therefore, all the dynamical participation ratios Y^^" 
are different from their static counterpart. The relative 
weights of the different visited sites are not given by the 
ratio of their Boltzmann weights. This result is impor- 
tant, since it was shown that, in the Sinai model, equilib- 
rium and dynamical participation ratios indeed coincide 
[45]. An interesting possibility, discussed in the context 
of glassy systems, would be that the Y^^'-'^ correspond to 
an equilibrium measure but at a different effective tem- 
perature. We shall discuss below that this is not the case 
either. 

Note that for walks in higher dimensions, d > 2, one 
can show rigorously that y^^" = Q [31], whereas Y^'^ 
are still given by Eq. (15)."'^ However, in this case, each 
site is visited a finite number of times, and therefore it 
could have been expected that the idea of partial equilib- 
rium over the set of visited sites would be quantitatively 
incorrect (although it is able to reproduce, at least quali- 
tatively, some non trivial dynamical correlation functions 
- [29] and see below). Another solvable case is the one di- 
mensional directed walk, where each site is visited once. 
In this case, ¥2^" can be computed exactly [33] and is 
found to be close to, but different from, the equilibrium 
value 1 — fi. 

The surprising aspect of our result in one dimension is 
that each site is visited, asymptotically, an infinite num- 
ber of times - a feature that, at least naively, should lead 
to partial equilibration. 

A different, but related, issue concerns the fraction 
fi(t) of the total time t a given particle has spent in 
the trap, and study the participation ratios of this 
quantity. In this case, we have found numerically that 
the different Yk are given by the equilibrium formula, 
Eq. (15). Therefore, we are in a situation where ergodic- 
ity is (weakly) broken: the relative time a given particle 



^The case d = 2 is marginal, but one still finds that Y^'"" = 
in that case [31]. 



6 



spends on different sites has not the same statistics as 
the relative fraction of particles that are found on the 
different sites at a given instant of time. Such a dif- 
ference between individual and ensemble measurements 
have been emphasized in a different context in [22], and 
recently observed experimentally [24]. 



C. Analytical calculation of the participation ratio 

Using the same procedure as for {p{x,t))r, one could 
try to compute Yfe {€) which is given by: 



/CO 
dx{p{x,tf), 
-oo 



(17) 



However, this calciilation reveals to be miich harder than 
for the case of {p{x,t))r. In Appendix C, we report a 
simplified calculation in the case k — 2, which aim is to 
argue that Y^^" is different from in the low temperature 
phase n < I, whereas it vanishes for /i > 1 . Although this 
last result has been rigorously proven in [30] , we want to 
introduce here a general method that could in principle 
yield bounds and approximations for Yk for any /x, and 
not only for < 1. We obtain the behaviour of Y^''^'" for 
/X ^ 1~ and find how Y2{t) vanishes as a function of t 
for /i > 1. To do this, we introduce a function R{t,t') 
through: 

/oo 
dx{p{x,t)p{x,t'))r (18) 
-oo 

(so that Y2{t) = R{t,t)), as well as its Laplace transform 

R{s,s'y. 



Ris,s') 



dt / dt'e-'*-''*'R{t,t') (19) 

"'0 



Using rather crude approximations, we obtain that, in 
the particular case where s = s' and /x < 1: 



R{s,s) ~ ^ 







(20) 



with a finite coefficient Rq. In order to interpret this re- 
sult, we assume that R{t, t') obeys, for large t, t', a scaling 
relation of the form: 



R{t,t') = Yi^"n(^^ 



(21) 



which we have confirmed using numerical simulations. 
Then one gets: 



R{s,s) - 
or, using Eq. (20), 



2K 



dyn 



(i+») 



2' 



(22) 



dyn 



Ro 



(23) 



Since the integral appearing in the above equation is con- 
vergent (because TZ{u) < 1), this result suggests that 
Y2^" is finite when /j, < 1. As furthermore we find that 
Ro vanishes linearly when /x — > 1~, we conjecture that: 



Y, 



dyn 



OC 1 — fi {fi ^ 1) 



(24) 



which is compatible with the numerics and also compa- 
rable with the equilibrium behaviour. The same level of 
approximation on Y^ also leads to a finite limit Y^^^, and 
to a linear temperature behaviour Y^^" oc 1 — /x (/(x ^ 1), 
so that one can reasonably guess that this linear depen- 
dence is valid for all fc > 1. 

The case n > 1 has also been studied; we find for 
1 < /U < 2 the new predictions: 



R{s, s) 



1 



5 — /J. 
S 2 



(1< M < 2) 
1 



(25) 



i?(s,s)~— (/x>2) 

S2 



which predicts that Y2{t) tends to zero as f^i"'')/^ when 
1 < /i < 2, and as fr^^^ whenever fi > 2. The last result 
is indeed expected: when the second moment of tp{T) ex- 
ists, diffusion is normal with no anomalous corrections. 
The probability that two particles starting at time at 
the same site happen to be again on the same site at 
time t decays as: ^{t) x l/(,{tY oc l/\/t. The result for 
1 < /i < 2 has been checked numerically for /x = 1.3, 1.5 
and 1.7 (see Fig. 6). 



Y,(t) 




FIG. 6. Plot ofY2{t) for /x = 1.3, 1.5 and 1.7 (Ml lines, 
from top to bottom), showing a power law decay compatible 
with the predicted behaviour t^^~^^^^ (dashed lines). Note 
that the sub-leading corrections become stronger when /x is 
close either to n = 1 or to fi = 2. 



7 



D. Partial equilibrium in a finite region 

We have seen that the dynamical participation ratio 
never reaches the static equilibrium value. Can one how- 
ever isolate a region of space, of size £{t) possibly much 
smaller than ^{t), such that inside that region equilib- 
rium is reached ? In order to test this idea, one can 
define a spatially restricted participation ratio Yk{£, t) in 
the following way as: 



(26) 



\i\<e 



function of £/£^{t). The collapse is rather good, showing 
that the size up to which the system is equilibrated grows 
as ^{t), which is thus the only dynamical length scale of 
this model. That Y2{£, t) is equal to the equilibrium value 
for £ < 4>S,{t), where </> is a small number, means that only 
the 'contemporary' processes concerning the largest scale 
^(f) are out of equilibrium. In a sense, this could have 
been expected. However, let us emphasize again that a 
simple description such as Eq. (11), which describes the 
lack of equilibrium on the scale of ^ through the form fac- 
tors gi{t), cannot explain the observed difference between 
Y^" and f/^". 



where Pi{t) is the probability that the walk is on site 
i conditioned to the fact that it is within the interval 



(27) 



Fig. 7 shows the numerical results for the following 
rescaled quantity: 



dyn 



AY2{£,t) 



Y2{£,t)-Y 



(28) 



such that ~ 1 for an equilibrated region, and by 
construction for £ :» ^ [t) . 




FIG. 7. Plot o{AY2{l,t) versus £ for n = 0.5 and t = 10^, 
lO"*, 10^ 10^ and 10^ so that ^{t) ranges from 10 to 215. 
Inset: the size £ of the window is rescaled by ^{t), and the 
resulting collapse of the curves shows that the equilibration 
length scale is of the order of a fraction of ^{t). The strong 
increase of AY2{£, t) for small £ is due to small size effects. 

The results are obtained with t ranging from lO'^ to 
10'', and 1^1 = 0.5. When t goes to oo at fixed £, Y2{i, t) is 
seen to converge to the corresponding equilibrium value 
(which depends slightly on £: Eq. 15 is only valid in the 
limit of large sizes). In the inset we show ^Y2{£,t) as a 



IV. GENERALIZED EQUILIBRIUM AND 
HALF-SPACE EXCURSIONS 



A. A functional relation between the Yit's 

We have seen that the dynamical participation ratios 
do not take their equilibrium value. Could one rede- 
fine an effective temperature /i such that all Yfe can be 
expressed as equilibrium values with this effective tem- 
perature? In order to test this idea, one can eliminate /x 
from the relation (15), and re-express all Yfe as a function 
of Y2. One finds: 



Yk = fk{Y2) = 



r(fc - 1 -f Y2) 

mnY2) 



(29) 



We have plotted in Fig. 8 Y.!^^"' versus Y^"'^ for several 
values of /U. It appears clearly that this relation is dif- 
ferent from the equilibrium one, shown for comparison. 
This rules out the possibility of defining a meaningful 
temperature from Y2^^ ■ 

The above relation between the Yk is known to be in- 
correct in other models, such as in the random map for 
example (see [41]). Inspired from the replica method, one 
can formally generalize Eq. (15) to: 



rdyn 



Yk^n — 



r(l - n) r(A; - At) 

r(fc - n) r(i - /x) ' 



(30) 



where n is the number of replicas (that must be set to 
n = to recover Eq. (15), see [41,44]). We can then ex- 
press Yk^n as a function of F2,„ for arbitrary n. We found 
that the value n = — 1 gives a reasonable account of the 
data for all ji values. Interestingly, this value n = —1 
was found to describe exactly the 'area preserving ran- 
dom map model' considered in [41] (where other models, 
corresponding to different negative values of n, where also 
studied). However, as we discuss now, a more precise in- 
vestigation of the problem shows that n = — 1 does not 
fully describe our results. 
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FIG. 8. Plot ofY,f''"{ii) versus Y2 ^"{fj,), parameterized by 
jj,. The equilibrium relation, corresponding to n = in the 
replica language, is shown for comparison (dotted line). A 
rather good account of the data is obtained using n = —1 
(dashed line). However, as we show below, a better descrip- 
tion is obtained using a mixture ofn = —2 cases, see Eq. (39) 
(full line). 

B. The dynamical distribution of weights 

Instead of studying all the different Yfe's, one can anal- 
yse directly the time evolution of the distribution of 
weights, (p{P,t), defined as (see also Eq. (12)): 

ip{p,t) = {Y,PHP-m)))r o<p<i (31) 

i 

For long times in an infinite system, this distribution is 

expected to reach a stationary distribution (pdyn 

(P). The 

inset of Fig. 9 shows <f{P, t) for three successive (large) 
times: t = 10^, 10* and 10^. All three curves collapse 
rather well, at least not to close to the 'edges' P = 
and P = 1, showing that we are close to the asymptotic 
distribution. (However, since the Yfc's are sensitive to 
the region around P = 1, the discrepancies at the edges 
explain why these moments converge more slowly). 

One can define a generalized distribution (fn,p.{P) as 
the one generating the Yfc,n(M)) which leads to the fol- 
lowing beta distribution: 

^"-'''^'° r(i-A)r(A-„) ^"''^-^''""" 

Following the same line of thought as in the previous sub- 
section, wc want to check more precisely if the data are 
compatible with n = —1. 




ii I \ I \ I \ I I I I 

0.2 0.4 0.6 0.8 1 



p 

FIG. 9. Fit of (p{P) with the n = -1 ansatz - Eq. 32 - 
using /i as a free parameter, for ji = 0.3 and n = 0.1. The 

Gt is worse and worse as /x is decreased. Note that for higher 
values of ji, the fits are better. Inset: tp(P,t) is plotted for 
different times t = 10^ (fuJJ line), 10"* (dashed line) and 
(dotted line) in order to evidence the convergence towards an 
out of equilibrium distribution ipdyn{P) (m = 0.5). 

Taking /i as a free parameter, one can try to fit the 
numerical data. For /i > 0.5 the fits obtained arc correct 
(data not shown). On the contrary, for /x < 0.5 the best 
fits appear to be quite unsatisfactory, in particular for 
/i = 0.1 (sec Fig. 9), showing that the n = — 1 ansatz 
does not fully account for the numerical data. This is 
due to the fact that (fin,ii{P) is a monotonous function 
whatever the value of fi, whereas ipdyn{P) becomes non 
monotonous for low values of fi. 



C. A simple analytical argument in the limit /u 

Although the out of equilibrium localization problem 
seems to be hard to tackle at finite temperature, a simple 
argument can be given in the limit /i ^ 0. This argu- 
ment accounts for the non trivial limits Y^^^'' and Fg'*^", 
which were found to be less than 1 for /x — > (see Fig. 5). 
If n is very small, then the largest trapping times accessi- 
ble after a given time t are strongly separated from each 
other. One can for example show that the distribution of 
the ratio R of the second largest time over the largest is 
p{R) = iJiR^~^, which tends to 6{R) when /z ^ 0. There- 
fore, in this limit, one can assume that the time elapsed 
before finding the deepest trap io occupied at time t is 
negligible compared to the time spent in zq.^ 



^The time spent on site io is actually much greater than the 
trapping time Tig , since the particle comes back to it a laxge 
number of times before finding a deeper trap (see section V). 
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So the problem beeomes equivalent to that of a random 
walk with no random potential, but with two absorbing 
boundaries (i.e. the traps with trapping time > t to the 
right and to the left of the initial site) at random po- 
sitions. If these absorbing sites are at distances, respec- 
tively Xr and xi from the initial position of the walk, then 
the probability to be absorbed by (say) the left boundary 
isp; = Xr/ (xi+Xr)- Since the initial site can be anywhere 
between these two sites with equal probability, one finds 
that pi is a random variable uniformly distributed over 
[0, 1]. Coming back to the trap model, it means that only 
two sites can be occupied, and the corresponding occupa- 
tion rates are uniform random variables. More precisely 
(p(P) can be written as: 

MP) =P t dpi [5{pi -P) + d{l -PI- P)] = 2P 
Jo 

(33) 



which leads for Y). to: 



no = ^^P-Vo(p)^p = ^ 



(34) 



For the particular cases fc = 2 and fc = 3 presented above, 
this gives Y2 = | and = 5, which agrees rather 
well with what can be extrapolated from the numerical 
data on Fig. 5. Moreover, Fig. 10 confirms that (p{P) 
converges towards ipo{P) when /x ^ 0, although rather 
slowly. Note also that ipo{P) can be written in the gen- 
eral form Eq. (32) given in the previous subsection, for 
the special choice n = — 2 and p, = —1. 



<p(P) 1 



|i=0.2 
H=0.1 
H=0.05 




FIG. 10. ip{P) for several small values of /u: n = 0.2 (full 

line), 0.1 (dot-dashed line) and 0.05 (dashed line). The curves 
seem to converge towards the asymptotic distribution ipo{P) 
(dotted line) for fi ^ 0, but the convergence looks quite slow. 
Note that P = and P = 1 are presumably singular points 
in this convergence process. 



D. Generalization of the argument to finite 
temperature 

The above argument can be reinterpreted in the fol- 
lowing way. In equilibrium, the zero temperature limit 
means that a single site dominates and contains all the 
probability weight. This is why Y^"^ 1 when /i ^ 0. 
On the other hand, in one dimension, the time needed to 
explore an interval of size L is 

which grows 

much faster than the time to exit the deepest traps 
(~ iyi/^") found in the interval. Therefore, if a deep trap 
is encountered in - say - the left region of the line, there is 
a siibstantial probability that the particle will not have 
time to explore the right region and equilibrate with a 
trap of comparable depth. This is the essence of the 
above argument: at zero temperature, the fraction pi of 
the weight captured by the left trap is uniform between 
zero and one, independently of the relative depth of the 
two traps. A simple way to generalize this argument 
for finite temperature is to assume that each half-space 
is independently equilibrated, and carries a total weight 
uniformly distributed between zero and one, as in the 
zero temperature limit. (Actually, as noticed above, we 
only need to assume that in each half-space the prob- 
ability distribution has the form given by Eq.(ll) with 
arbitrary factors gi: this does not affect the asymptotic 
shape of ip{P) = ifieq{P)-) Denoting by ip{P,pi) the dis- 
tribution restricted to the left half-space, normalized to 
Pi, one has for < P < pi: 



^{P,Pl) = 



Pi 



r(i-M)r(M) 



p-"" {pi - py-^ (35) 



Averaging over pi with a uniform weight, and taking into 
account the right half-space, leads to the following pre- 
diction for lyj(P): 

^{P) ~ <p*(P) = 2 / dpi v{P,pi) e{pi - P) (36) 
Jo 



^*{P) 



r(i-M)r(i + M) 



+ 



2/i 



r(i - Ai) r(2 + m) 



p^-^ii-p)" 
p-ft(i _ p)i+^t 



(37) 



Although we do not have any interpretation for this, one 
can notice that ip* can be written as a superposition of 
distributions <^-2,/i, with two different values of fl: 

if* (P) = (1 - m) </'-2,^-i (P) + M ¥'-2,^ (P) (38) 

with (.pn^fi defined in Eq. (32), and in agreement with 
what was found for /x ^ 0. This prediction is com- 
pared with the numerics in Fig. 11, for several values of 
(/i = 0.1, 0.3, 0.5). The agreement is rather good; note 
that no fitting parameter is used here. We note that the 
participation ratios Yk obtained using <^*(P) are simply 
proportional to the equilibrium values: 
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n* 



1 



(39) 



As shown in Fig. 5, this relation accounts quite well 
(but not exactly) for the data. However, the distribution 
W{pi) of the weight pi = 1 — Pr carried by one half-space 
is not found to be uniform as we assumed, except for 
/i — > - see the inset of Fig. 11. Surprisingly, if one redo 
the above computation with a humped shaped distribu- 
tion W{pi) = Alpi{l—pi)]'^ , the resulting ip*{P) does not 
fit the data as well as the above form, which corresponds 
to (T = 0. 



<P(P) 2 




FIG. 11. Comparison between ip{P) obtained by numeri- 
cal simulations (symbols) and f*{P) given by the argument 
developed in the text (lines). The agreement appears to be 
quite good, bearing in mind that no free parameter is used. 
Inset: distribution W(pi) of the probability weight pi carried 
by one half space, for different /i. This distribution appears 
to be non uniform (except for jj, — > 0), at variance with the 
hypothesis underlying Eq. (36). 



E. Discussion 

The physical picture emerging from the last sections 
is the following: for the purpose of understanding global 
localization quantities, one can reasonably consider that 
space is split into two half lines, and that each one be- 
haves as if it was independently equilibrated, but out of 
equilibrium with respect to the other. On the other hand, 
we have seen in section HID that length scales much 
smaller than ^ could be considered as equilibrated, and 
that departure from equilibrium comes from the largest 
length scale |a;| > ^. The equilibration of the small scales 
cannot be accounted for by the previous argument. All 
these observations suggest, in order to get a consistent 
picture, that space may be actually divided into three 
regions, an equilibrated domain centered on the origin 
and two quasi equilibrated regions on each side, and that 
each part of space is not equilibrated with the others. 



An artificial remedy to this lack of equilibration be- 
tween the different regions is to allow the particle to 
make long jumps. We have therefore added 'links' be- 
tween the sites x and —x, such that the probability to 
hop directly from x and —x decays as x"''. When p is 
large enough, the dynamical participation ratio Y2{t) is 
a decreasing function of time, and seems to converge to 
y^dj/n^ ^hen p < on the other hand, Y2{t) is seen 
to reach a minimum and to increase back towards the 
equilibrium value Y^'' . We have however not checked in 
details whether Y2{t) indeed converges towards 1^2'^' for 
all p < Pc, but only wanted to illustrate that the differ- 
ence between Y^^^ and Kg*^^ is due to the scarcity of the 
links between the different sites for the one dimensional 
lattice. 



V. CORRELATION FUNCTIONS, AGING AND 
SUB-AGING 

A. Motivation 

Let us now turn to different correlation functions that 
one can define in order to probe the peculiar aging prop- 
erties of this model. Since the largest encountered trap- 
ping time during tw scales as t!^ with z/ = < 1, 
one would naively expect that two time correlation func- 
tions vary on a time scale ^ f^. This would correspond 
to 'sub-aging' behaviour, where the effective relaxation 
time grows less rapidly than itself. 

This is indeed the case for the probability 11 (i^ +t, i^) 
of not having jumped at all between time and + t. 
This correlation function was computed numerically in 
[29], and was found to scale very accurately as: 11 (<^ -|- 
t,tw) = T^ii/tw)- The shape of the scaling function was 
compared to the prediction of an approximate calcula- 
tion where one assumes 'partial equilibrium', i.e. that 
the probability to find the particle in a trap of depth 
T after time is equal to the equilibrium probability 
within a region of size i,{tw). This approximation pre- 
dicts a power law behaviour for 7r(s) both for small and 
large s, with exponents that agree with their numerical 
determination. The detailed shape of 7r(s) however de- 
parts from the numerical results, which is expected. The 
success of the partial equilibrium assumption here is due 
to the fact that n(i^ -I- t, i^) only depends on the aver- 
age probability to occupy a site, and not on higher order 
correlations such as needed to compute the participation 
ratios Y^. 

Perhaps surprisingly, different correlation functions 
may exhibit a completely different aging behaviour. Con- 
sider the probability C{tyj + t,ty,) that the particle oc- 
cupies the same site at time +t and at time t. Obvi- 
ously, C{tw + t,tw) > n(i^ + t, t^). But in this case, it 
was shown rigorously in [30,31] that C{tyj +t,tw) scales 
as a function oit/tw, and not as t/f:^. This means that 



11 



even if the particle has almost certainly jumped away 
from its starting point after a time <C tw, it has re- 
turned there even after a time of order so as to make 
C{2t^,t^) = 0(1), whereas Il{2t^,t^) 0. This dilfcr- 
ence is not intuitive a priori, in particular because one 
knows that once the particle has left its initial trap af- 
ter a time ^ t^, it takes on average an infinite time to 
get back there since the walk is one dimensional. But if 
C{tw+t, tyj) is to decay on the scale ty,, it means that the 
probability not to find the particle on its starting point 
after a time t much greater than but much less than 
tyj must tend to zero when — > oo. The fact that the 
particle jumps back and forth a large number of times 
between t!^ and could thus a priori lead to an inter- 
esting behaviour of C{tyj+t, tyj) in the short time regime 
t/tw <C 1 (which was not investigated in [30]). For ex- 
ample, one could find, as in [29], different 'time domains' 
t ^ f^, t ^ t^, etc., where the correlation function has 
a different analj^tic behaviour. This is the issue that we 
discuss below. 



B. Analytical arguments 

The difference of scaling between 11 and C can be qual- 
itatively understood as follows: for C(tu> + t,tw) to de- 
cay to zero, one has to wait until the region probed by 
the particle at time + t is much larger than the ini- 
tial region where it was located, i.e. a time t such that 

^{tw + i) > ^{tw)- But since ^(i^) ^ tiu"" , the time 
needed for C to decay to zero is necessarily of order ty,. 
(Note that this argument does not hold for n(t^ t^,), 
which only requires the particle to hop once out of its 
initial trap.) 

In the same spirit as Ref. [29] for Ii{ty, + t,tu,), but 
using a slightly difi'erent method, one can try to give an 
approximate calculation of C{tyj+t,tyj). The first step is 
to introduce the dynamical distribution of trapping times 
p{T,tw), assumed to behave as: 



(40) 



This encodes the fact that typical trapping times arc of 
order f^. If one assumes that short time scales (t <C tj^) 
are equilibrated, whereas large ones {t^ V^) are still dis- 
tributed according to the a priori distribution (this can 
be rigorously proved in the fully connected trap model), 
one obtains the following asymptotic behaviour for (j){z): 



4>{z) ^ ^ 







(j){z) ~ 



Using the relation: 



7o 



-|-oo 



(41) 
(42) 



/OO 
dTp{T,ty,)e-*/'' (43) 

one can easily deduce from Eq. (41) the short and late 
time behaviour of Il{tyj + t, tyj): 

n{tyj + t,tyj)c^i- T^r(M) - i « C (44) 
n(t^ + t,i^) ~7oor(At) (tit 



i»C (45) 



in agreement with the results of Ref. [29] , and with the 

numerics (sec below). 

Turning now to C{tyj + 1, iiu), one has to take into ac- 
count the fact that when a particle leaves its trap, it will 
come back a large number of times before really escaping. 
We thus propose the following approximation: a particle 
will be considered to have truly left its initial trap if it has 
encoimtered a deeper trap during its excursion out of the 
original trap. Given a trapping time r, the probability 
that r' > T is given by: 



P^r'>r)=r^ = ^ 



(46) 



So the probability p{i, t) that the first trap encountered 
with a trapping time larger than r is found at a distance 
£ is: 



p{i,r) = P{r' > T) [P{r' < T)t' = (l - ^ 



For large r's, one has: 

p{e,T) 



(47) 



(48) 



Note that we only give here a scaling argmnent, and that 
corrections coming from the fact that there is a deeper 
trap on both sides are neglected. Conditioned to the fact 
that the deeper trap is situated at a distance the par- 
ticle has a probability j to reach it, once it has jumped 
out of its initial trap. Therefore, the escape rate can be 
written as: 



1 



(49) 



The correlation function C(t^ + t^tyj) is then given by: 

dTp{T,tyj) j d£p(£,T)e-™(^'^)* 

(50) 

After a few changes of variables, and using the scaling 
relations Eq. (41), one finds the following short time and 
late time behaviour for C{tyj + t,tw): 
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City) H~ ^) ^iu) — Q 



^71 



where the constants Cg and c; are given by; 

2 



Ci = 



1-M 



1 



1 + M 



(51) 
(52) 

(53) 
(54) 



These values for the short time and late time singular- 
ity exponents, have, to our knowledge, not been reported 
before, although they should in principle be contained in 
the analysis of [30] . We now turn to a numerical investi- 
gation of these asymptotic predictions. 



C. Numerical results and multiple time scales 

1 

Fig. 12 displays l — Yi{ty,+t, t^) as a function of i/ti"^" , 
and C{tw +t,t,u]) as a function of for different wait- 
ing times {tyj = 10"^, 10^, 10^ and 10^) and at tempera- 
ture M = 5- The collapse is very satisfactory, confirming 
the validity of the predicted scaling relations. 




FIG. 12. Plot of 1 — Il{tw + t,tw) (increasing curve) ver- 
sus t/t'^ and C{tw + t,tw) (decreasing curve) versus t/tw, for 
fj, = ^, and V = = 2, Tie scaling relations arc very well 
satisfied, at least in this time window. Symbols refer to the 
same waiting times for the two curves: tw = 10^ (+), 10* (o), 
10^ (o) and 10® (>). 



On the other hand, a similar plot of ln(l — C) as a 
function of hi{t/tw) is less convincing, which could be 
the sign of multiple time regimes (as was the case in [29] , 
where similar 'non scaling' features actually suggested 
such regimes). A way to investigate this issue is to study 
the function g{a,tyj) defined as: 



g{a,tw) = - 



in[i-r(/„ + /;^/,„)] 



(55) 



If this function has a limit 500(0;) when t^} — > 00, it means 
that in the time domain where t ^ f^, the probability 
1 — C that the particle has escaped from its starting site 
decays as tw^°°^"^ for large tw Prom the t/tw regime 
established by [30], wo already know that .9oc(l) = 0. 
If 1 — C{tyj + t,tw) behaves as (t/tw)^ even for i ~ 
with a < 1, then one should observe goo{oi) = A(l — a). 
Any departure from a linear function goc (a) would signal 
multiple time regimes; in particular, for the model con- 
sidered in [29] in = 1 where two sub-aging exponents 
z^2 < t^i < 1 appear, one finds that the function (700(a) is 
piecewise linear in the intervals [0,1^2] and [v2,vi], with 
different slopes. One also finds that goo{'^2) — 9oo{^t)^ 
and goo{ct > i^i) — 0. In this case, the change of 
slope indicates the presence of a characteristic time scale. 
One could imagine more complicated 'multi-scaling' sit- 
uations where goo (a) is a non trivial curve. 

We have first tested this procedure on U{tiu + t,tw), 
defining in the same way a function g*{a,tw) associ- 
ated to n. In this case, a single sub-aging scaling is 
expected, with v = and A = 1 — /i, which leads to 

5oo(a) = (1 — m)(^ ~ o). The function g*{a,tw) is plot- 
ted for different values of t^ (namely t^ = 10'', 10^, 10^ 
and 10^) in Fig. 13, for /i = i. On general grounds, one 
expects finite time corrections to g'^ia) that decays as 
l/\ntw and l/t^,- Using these corrections, one can very 
satisfactorily extrapolate g^{a,tw) to a function which 
is very close to the expected result g'^ict) = 5(| ~ Q^) 
(see Fig. 13). To be more specific, we used the following 
functional form for the extrapolation: 



g*{a,tw) = g*^{a) 



hitu 



c{a)t 



-7(a) 



(56) 



Let us analyse in more details the short time behaviour 
of these correlation functions. When plotting ln(l — 11) 
as a function of ln(i/t5^), the scaling is still quite well 
obeyed and in good agreement with the theoretical pre- 
diction 1 — n ~ {tltwY~^, obtained in [29], up to small 
time corrections that vanish only when T^t ^ 1. 



where 3^(0;), c{a) and 7(a) are fitted for each value of 
a, and 6 is a fitting coefficient independent of a, since the 
l/\ntw correction is expected to come from the prefactor 
of {t/t'^y~'^ in the short time expansion of the correla- 
tion function. Therefore the value of b was fixed from the 
direct power law fit of the short time regime of 1 — 11. 
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a 



FIG. 13. Function g*{a,tiu), associated to n(fu, + t,<u,), 
with M = I and ty, = lO", 10'''' , lO'' and lO'^. The function 
g^{a) (full line) is expected to be g^{a) = i(| — q) (see 
text). The infinite time extrapolation g*{a,tiu) (o) agrees 
very well with the prediction, with small discrepancies near 
a = and a = | , where finite time effects are stronger. 

One can now apply the same procedure to C{tw+t, tw). 
The results are shown in Fig. 14, using the same con- 
vention as for Fig. 13; g{a,tu,) is represented for the 
same waiting times as g*{a,tw). Interestingly, although 
finite time corrections are strong, the extrapolated re- 
sults are in good agreement with our analytical predic- 
tion Qooict) = A(l — a), with A = (1 — /i)/(l + n) = 1/3, 
at least when a G [0.2, 0.8]. This suggests that a unique 
time regime t ^ tw is relevant for C{tw +t,tw), although 
we know that the time scale <C tyj governs the evolu- 
tion of n. Note that the l/\ntw corrections are weaker 
than in the previous case, and one is almost dominated 
by power law corrections. This is due to the fact that the 
prefactor of (t/tw)^ in the short time regime happens to 
be close to 1 here (and hence the parameter h is small), 
whereas the prefactor of {t/f^)^ was about 0.57 for 11. 

So, what happens to the particles that have left the 
initial trap after a short time and took a very long 
time to come back? The probability that a particle 
leaves the trap exactly at t' is d[l — U{tm,tw + t')]/dt' ~ 
t'~'^ /t^^~'^\ If the sample was not disordered, the prob- 
ability that it has not returned to the origin after time 
t — t' decays as (t — Because of the long trap- 

ping times, this probability actually decays slower, as 
{t - These particles contribute to 1 - C, as: 

Jo tT ^' 

Choosing t = t^, we find that the contribution of these 

'early birds' to 1 — C is a factor tw°'^ /(i+a*) gjj^aller than 
the contribution computed above, Eq. (51), and are thus 
negligible in the large tw limit. 




a 



FIG. 14. Function g{a,tw), associated to C{tw +t, tuj), for 
/i = 5 and with the same waiting times as for g*{a,tw). 
The argument developed in the text - Eq. (51) - predicts 
goo{cx) = |(1 — a) (full line). Although finite time corrections 
are stronger than in the previous case, the infinite time ex- 
trapolation gexia) (<>) agrees well with the prediction, at least 
for 0.2 <a< 0.8. 

This simple estimate shows (i) why finite time cor- 
rections become large when a — > 0, (ii) that power-law 
corrections to g(a,tw), as the one used for our extrapo- 
lation, are indeed expected and finally (iii) justifies why 
the time scale does not appear in C{tyj +t,tyj). 




FIG. 15. Plot of the late time behaviour of'n(tw+t, t^) ver- 
sus t/t™, and C(tw -\- t,tw) versus t/tw, for tn, — 10^ and 10^ 
respectively, and fj, = 0.5. Both correlation functions exhibit 
a power law behaviour at large time, and the exponents agree 
well with the predicted values — i and — ^ (see Eqs. (44) and 
(51)). The corresponding slopes are shown in dotted lines, as 
a guide to the eye. 

The late time behaviour of the two correlation func- 
tions can also be tested numerically. Fig. 15 shows 
Il{tw + t, tw) versus t/f^, and C{tw + t, tw) versus t/tw^ 
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for = 10^ and 10^ respectively, and ^ = 0.5. The two 
correlation functions behave as power laws at large time 
t, and the exponents are in good agreement with those 
predicted by Eq. (44) and (51), shown for comparison 
(with a arbitrary prefactor), at least for this particular 
value of /U. Note that since we do not know the constants 
7o and 700, we cannot test the values predicted for the 
prefactors. 



search, and that some of our results, in particular con- 
cerning asymptotic estimates, can be proven to be exact. 

ACKNOWLEDGMENTS 

Fruitful discussions on this model with G. Ben Arous, 
J.-M. Luck, C. Newman, P. Maass, M. Mezard, C. Mon- 
thus and G. de Smedt are gratefully acknowledged. 



VI. CONCLUSION 

In this paper, we have studied in details the one di- 
mensional exponential trap model, which exhibit a phase 
transition between a high temperature diffusive phase 
and a low temperature sub-diffusive phase. We have ob- 
tained numerically and analytically the shape of the av- 
erage diffusion front in the sub-diffusive phase. Although 
based on an approximation valid only near the dynamical 
transition, our calculation provides several predictions on 
the asymptotic shape of {p{x, t)) which are in excellent 
agreement with the numerics. It would be interesting to 
see whether these predictions are actually exact. 

The central result of this study concerns the localiza- 
tion properties. We have found that the dynamical par- 
ticipation ratios are all finite, but different from their 
equilibrium counterparts, even allowing for the existence 
of an effective, dynamical temperature. This is surprising 
because since each site is visited a very large number of 
times by the random walk, one could have expected that 
a partial equilibrium sets in within the limited region of 
space explored by the walk. Our detailed study of the 
distribution of dynamical weights shows that this is not 
the case. We have argued that this can be interpreted in 
terms of an effective 'fragmentation' of space in two half 
lines (or even three domains), with a restricted equilib- 
rium within each region, independently of the others. 

Finally, we have stuclic;d two different two-time corre- 
lation functions, which exhibit different aging properties: 
one - Il{tyj+t, tyj) - is 'sub-aging' whereas the other one - 
(^{tw +i) tw) ~ shows 'full aging'. We have given intuitive 
arguments and simple analytical approximations that ac- 
count for these differences. We have obtained new predic- 
tions for the asymptotic (short time and long time) be- 
haviour of the scaling function associated to C{tyj +t, tw), 
which are found to be in excellent agreement with the nu- 
merics. Since two time scales {t'^ and t^) appear in this 
model, one can wonder whether the short time behaviour 
of C{tw + t,tw) exhibits a non-trivial, multiple time scal- 
ing. A careful numerical investigation of this issue leads 
to a negative answer, although strong finite time correc- 
tions are expected. 

Since this one dimensional model is currently of inter- 
est to the mathematical community, we hope that the 
study presented here will motivate further rigorous re- 



APPENDIX A: THE AVERAGE DIFFUSION 
FRONT 

Formulation of the problem 

The explicit calculation of {p{x,t))r is reported is this 
appendix. Two different averages will be introduced, the 
average over the random walks (...)„, and the average 
over the quenched trapping times (. . .),-. We consider 
here a slightly modified version of the model, in which the 
particle stays on a site a time exactly equal to t{x) (con- 
tinuous notations are used, so as to facilitate the contin- 
uous space limit) rather than exponentially distributed 
around this value. We expect, and have checked numer- 
ically that this is irrelevant for the shape of {p{x, t))r at 
long times. 

For a given sample of the disorder (quenched trapping 
times), we can decompose the probability p{x,t) for the 
walker to be on site x at time t into a sum over the num- 
ber n of steps: 

00 00 

*) = X] ^ '^{^Xn,xl{tn < t < tn+l))w 

71=0 71=0 

(58) 

where /(t„ < t < <„+i) is the characteristic function of 
the interval [tn, tn+i], equal to 1 of t belongs to this inter- 
val, and otherwise. In order to simplify the notations, 
we introduce In{t) = I{tn < t < tn+i)- Now, averaging 
over the disorder: 

00 

{p{x, t))r = ^((<5x„,x^(t)).^)x (59) 

The key-point is that we can permute the two averages, 
and perform first the average over the disorder for a given 
walk. Introducing the average (. . .)n,x over the n steps 
walks ending on site x, we get: 

00 

(p(a;, t))r = Q{x\n){{In{t))r)n,x (60) 

where q{x\n) is the standard probability for the random 
walk to be on site x after n steps: 

q{x\n) = -L= e-^'/2" (61) 
v27rn 
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for large n. Taking the temporal Laplace transform C of 

{p{x,t))r: 



{p{x,s))r = '^q{x\n){{in{s))r) 
n=0 

which requires the calculation of In{s)- 



(62) 



Inis) = I e 



-'^*dt = - e-'*- [1 - e-«(*''+i-*"->] (63) 



since tn+i — tn = t{x). and t{x) is at most of order s^'"' 
when s — » 0, so that st{x) should be small. For a given 
walk W ending on site x after n steps, the time t„ can be 
decomposed into a sum over the different visited sites: 



tn = J2N'w{x')T{x') 



(64) 



where Afwix') is the number of visits of the site x' by the 
walk W. Now /„(s) can be averaged over the disorder: 



(/n(s)). = (Te-«^-(-)-). Y[{e 



)r (65) 



Averages of the form (e °''^)r or (re '*^)r are easily cal- 
culated, in the limit a ^ 0: 



l-ca''~e- 



d 



da 



(66) 
(67) 



with c = r(l — fx), so that (J„(s))r reads: 

{in{s))r = CM[sM^(a:)]''-'e-'^«''[S»'^'^(^')1 (68) 

Next, one has to average over all the walks W ending on 
site a; in n steps. Since this part is the hardest one of the 
calculation, one has to resort to a simple approximation 
scheme, valid in the vicinity of some specific value of fj,, 
namely jj. close to 1 in the following. 

An approximation for /u — > 1 

A simple approximation consists in performing the av- 
erage (. . .)n,x of the right hand side of Eq. (68) by re- 
placing A/V (2;') with A4(;c', n) = {f^wix'))n,x' 



{{In{s) )t )n,x = {Cfj. [sAfwix)] 



x-r]\ 

/n,x 

~ cfx[sNx{x,n)f-^ e-=^''[S^'^==(^''")1 (69) 

This approximation is expected to be correct, at large 
times, for close to 1 (/z < 1), since it is exact for /x = 1. 
Note however that for /z = 1, Eq. (66) is no longer valid. 



and logarithmic corrections come into play. Turning to 
J\fx{x',n), it can be shown that for large n, a scaling re- 
lation holds: 



Afx{x',n) = sJ^F ( 



with F{y, z) given by the following integral: 

du 



F{v,z) 



\/2^ Jo 



. g 2 sin nu 



(70) 



(71) 



3m TTU 

The different factors in Eq. (69) can then be evaluated: 
f +00 



/i-oo 
-00 



(72) 



X 



introducing G{z) = F{v, z)^dv. 



X X 



M-l 



(73) 



with Fl(y) = F(v,v). From Eq. (62), a continuous space 
limit can be obtained, introducing a continuous scaling 
variable A through the natural scaling relation n = Xx^: 



{p{x, S)}r = x"^ 



roc 


■ g-l/2A - 


/ d\ 




10 


\x\V2^_ 



(74) 



X e 



where the parity of G{z) and H{z) has been used. Group- 
ing together the factors, and introducing the scaling vari- 
able 77 = |a;|(-^+'^)/^s, we get: 



roo 

{p(x,s))r = \x\^^^-^V^-^ dXX^- 

Jo 



(75) 



Note that this expression is compatible with the expected 
scaling form: 



{p{x,t))^ 



(76) 



where ^ is the dynamical length scale appearing in the 
model: ^ ~ One can indeed rewrite the previ- 

ous relation in the following way: 
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{p{x,t))r = ^,g ^ 



(77) 



Taking the Laplace transform with respect to t yields: 

{p{x,s))r = \x\''^g{v,) (78) 

where g is the Laplace transform of g, which is of the 
form Eq. (75). So one deduces that the scaling function 
g{v) is given by: 



9 iv) 



IJ,C 



Jo 



(79) 



._i__c<a(i+'''/=^g( 



•/J J 



Asymptotic behaviour of the scaling function f 

Now we focus on the asymptotic behaviour of g{r]) for 
large r], which gives the spatial tails of the distribution 
{p{x^t))r- When 7] oo, the above integral is domi- 
nated by the small A region, which means that one needs 
to know the asymptotic large z behaviour of H{z) and 
G{z). After a few lines of computations, we finally find: 

H{z) ~ - G{z) ~ z^-f" z^oo (80) 
z 

The large r] behaviour of g{r]) can be then obtained from: 



lie 



r+co 

rj"-^ / dAA^-ie-^-^""^". (81) 
Jo 



The inverse Laplace transform can be computed using a 
saddle-point method. One finally finds for the large \x\ 
behaviour (or large \(\, with ( = xt~T^): 



with: 



fo 



/(C)«/oo|CI 



/^r(i - m) 



(82) 



b = 2-^r(l - /x). (83) 



One can also look at the limit (, ^ 0. Starting from 
Eq. (79), we have to calculate the small z behaviour of 
G{z) and H{z), which is simple here since these func- 
tion have a finite limit in 0, denoted by 5o(m) ho 
respectively: 



I-oo ( \/27r Jo 



du 



2tt Jq \J sin TTW 
du ^ ' ' ^ 



r(ir 



\/sin TTU 

So we have to compute the following integral: 



(84) 
(85) 



gin) = /ir^^-i / d\ x^-h-^-osoWV^^'^''''' 

V 27r Jo 

(86) 

For ?7 0, this integral is dominated by the large A 
region. We finally find: 



g{ri) = Ar] ^+f^ rj 



(87) 



with: 



1 2/x 



-^r(l - „)rhg,(^)-T^hr^T (-^ 



(88) 



Taking the inverse Laplace transform, this term gives the 
value of {p{x = 0, t))r which is proportional to . In 

order to get the spatial dependence, the next term of the 
expansion must be computed. After a few changes of 
variables and asymptotic estimates, we get: 



/(c) = /o-/iicr c-o 



with: 



(l + A<)V27r 



r(i^/i)^r(^^ 



(89) 



(90) 



(91) 



The next sub-leading term can also be computed, and 
is found to be of order C,^ for > 1/2, and (^^^^f^ for 
M < 1/2. 



APPENDIX B 

We give here some technical details about the numer- 
ical simulations. A number of independent 'walk- 
ers' (or 'particles') are simulated one by one, for a given 
sample of the quenched energies {Ei}. A walk is sim- 
ulated as follows: the trapping time on site i is cho- 
sen randomly from an exponential distribution of mean 
Ti = exp{Ei/T), and the walker then chooses at random 
between the two neighbouring sites, with equal probabil- 
ity. Then the desired quantity is computed for this partic- 
ular sample, and eventually averaged over a number Ng 
of samples. Moreover, in order to facilitate comparisons 
between different runs, we took each time the same disor- 
der samples, by choosing the same set of 'seed' numbers. 
For out of equilibrium simulations, where the number of 
sites L = 2N + 1 is supposed to be infinite, we used peri- 
odic boundary conditions, with usually N = 10^ except 
when long times were required, in which case N = 10* 
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was instead considered (for instance in the computation 

of nw)- 

Error bars are estimated by running several simula- 
tions with the same Ny^ and Ns, varying only the 'seed' 
numbers. The fluctuations between the different runs 
leads to an estimate of the standard deviation. The num- 
bers 7Vu> and A^^ were chosen so as to get small enough 
error bars, as far as possible, taking into account the time 
t that we need to reach as well as the computational time. 
It should be emphasized that the amplitude of the fluctu- 
ations depends a lot on the computed quantity. Correla- 
tions functions are easy to compute, and Ny, = Ng = 10^ 
is enough to get a relative standard deviation which is 
less than 10~^. Quantities related to localization prop- 
erties fluctuate more; for integrated quantities like the 
participation ratios, we took N^, = 10^ and Ng = 10'* 
(except for small size systems where Nyj = Ng = 2.10^ 
was used instead), which was a good compromise in order 
to obtain both a good standard deviation (less than 10~^) 
and long enough times (for instance t = 10^ for fj, = 0.5). 
For distributions like (p{P) and {p{x,t)), one needs to 
average over a larger number of samples, in order to get 
smooth enough curves (with fluctuations between differ- 
ent runs less than 5.10"^). So ip{P) was simulated using 
= 10^ and Ng = 10^, whereas {p{x, t)) was computed 
with iV^ = 10^ and Ng = 10^ Note finally that Y2{£,t) 
was simulated using Ny, = 10^, and Ng = 10^ (except 
for t = 10^ and 10^, where A^,„ = 5.10* and 10^ respec- 
tively), since we need a good statistics on the sites with 
small trapping times so as to avoid large fluctuations at 
small i. 



APPENDIX C: CALCULATION OF THE 
PARTICIPATION RATIOS 



The analytical calculation of Y2{t) appears not to be 
easily tractable, and the aim of this appendix is to ar- 
gue for the existence of a flnite limit y^**^" of Y2 (t) when 
f — > 00 for /i < 1, and try to extract some information 
on the behaviour of y^'*^" as a function of fi, in particular 
for /X close to 1. 

The participation ratio Yzit) is given by the integral 
{p{x , t)"^) r ■ The quantity {p{x,t)'^)r can be com- 
puted following the same lines as for {p{x, t))r- It will be 
useful to introduce a two-time quantity Q{x, t, t') defined 
as: 



Q{x,t,t') = {p{x,t)p{x,t'))r 



(92) 



Defining R{t,t') = J^^dx Q{x,t,t'), one has Y2{t) = 
R{t,t), and turning to the Laplace transform: 

/•OO /'OO 

R{s,s')= (it dt'e-'*-''*' R{t,t') (93) 
Jo Jo 



Now a reasonable assumption (that has been checked nu- 
merically) is that for large t and t' , t > t' , R{t, t') becomes 
a function of 4: 



R{t,t!) = Y^y''n(^-}j 



(94) 



This follows from the similar behaviour of the correlation 
bmction C{t,t') studied in section V. It is interesting to 
restrict to the particular case s = s': 

dte-'^'+''^n(j\ (95) 



^(s,s) = 2y/^" r dt' f 
Jo Jt 

= 2 Y^y^" j dt' j dut'e-^*'(i+")7e( 



or: 



R{s,s) 



2Y. 



dyn poo 



du 



n{u) 



(96) 



So a finite limit for Y^it) corresponds to R{s,s) ~ s~^, 
and this is what we shall try to show in the follow- 
ing. Coming back to Q{x,t,t') and decomposing over 
the number of steps, one has: 



Q{x,t,t') = ^p{x,n,t)p{x,n',t') 



(97) 



Averaging over the disorder yields: 

{Q{x,t,t'))r = Y,q{x\n)q{x\n') (98) 

X (7(i„ < t < i„+i)/(C <tf< C+i)) T,{n,x),{n' ,x) 

For given walks W and W, and a given sequence of Tj, 
let us introduce Kn,n'{s, s') defined by: 

poo poo 

Kn.n'{s,s')^ dt dt'l{tn<t<tn+l) (99) 

Jo Jo 

X /(C <t'< C+i) e-^*-^'*' r{xfe-^'--^'''^' 

assuming again that t^+i ~ t^ — t'^,^^ — t'^, = t{x) 
(i.e. trapping times are fixed rather than exponentially 
distributed) and that st{x) and s't{x) are both much 
smaller than 1 , which means that the maximum trapping 
time encountered is much smaller than the times t and t' 
considered. Introducing the following decomposition: 



tn = ^Afw{y,n)T{y) 

y 

y 

Kn,n'{s,s') reads: 



Kn,n'{s,s') = t{x 



:2 —[sAfw (x .n)+s' N'nri (x,n')]T(x) 



(100) 
(101) 

(102) 



xne 

yjtx 



-[sN'w(y,n)+a'Nwi(y,n')\T(y) 
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Averaging Kn,n'{s, s') over the disorder, one has: 

{Kn,n'{s,s'))r=tiT{2-fJ.) (103) 

X [sJ^wix,n) + s'J^w'{x,n')]^~^ 

X Q-c[sAfw{v,n)+s'J\f{y{y,n')Y 

One can now write R{s, s) as: 
R{s, s) = 2iiV{2 - //) J2 q{x\n)q{x\n')s^'-'^ (104) 

X n<n' 

X [Mw{x,n) + Uw'{x,n)]^~'^ 

X f,-cs''T.y[j^w{v,n)+Myy,{v,n')Y 

We now turn to continuous hmit, and replace as in Ap- 
pendix A AfwiVi by its average value y/nF (^-^, ■ 
At this stage we drop order unity constants since we shall 
make several rather crude approximations in the follow- 
ing. Introducing the new variable /? through n' = f3n, 
one gets: 

R{s,s) ~ s'^-^ r dx r dn r ^e-^^^+^^ (105) 

J-oo Jl J I VP 

/f3nH ' 



X e 



So this simplified calculation is consistent with a finite 
Y^ "'^' given by: 



Y^y- oc — ^ r § r dx 

2 r(i-/x)A (i^J_^ 



TT I X 



^-2 



(109) 



where we have used c = T{1 — fi). 

Since lim2^oG'(2;) = goifJ^) 1 when n ^ 1 and 
lim2^o-ff(-2) = ho is independent of fi, the only strong 
dependence upon ^ comes from r(l — /j,), which suggests 
that: 



y/^"~(i-M) 



(110) 



This results is compatible with numerical data, and is 
comparable to the corresponding result in equilibrium 
(in this case, Fj^^ = 1 — for all /x < 1). For 1 < /x < 2, 
one can easily show that, under the same assumptions, 
Y2{t) ~ l/t('^-i)/2 when t -» oo, whereas Y2{t) ~ 1/Vi 
for /X > 2, which is the expected result. Turning to Ys, 
wc have checked that the same calculation also leads to 
a finite limit, and to a linear behaviour with respect to 
/U when ^ 1. 



Now, because of the factor the integral over /? is 
dominated by the large /3 behaviour, which means that 
as a first step one can neglect terms like F (^-^, 

compared to y/PF (^^^^i • Rescaling also x and y 
using the new variables x = and y = -^j^ yields: 



R{s,s) 



X e 



oo POO 

"-^ I dx J dfi(3' 



dnn^ (106) 



VP) 



One can change variable in the last integral over n, let- 
ting: 



1 



n = 



2m 



(107) 



Then the integral over v becomes Jq°° dve = 1 (for 
s ^ 0), and R{s, s) reduces to: 



H-2 



(108) 
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